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Abstract. We consider the transport properties of topological insulators surface 
states in the presence of uncorrelated point-like disorder, both in the classical and 
quantum regimes. The transport properties of those two-dimensional surface states 
depend strongly on the amplitude of the hexagonal warping of their Fermi surface. It 
is shown that a perturbative analysis of the warping fails to describe the transport 
in experimentally available topological insulators, such as Bi2Se3 and Bi2Tc3. Hence 
we develop a fully non-perturbative description of these effects. In particular, we find 
that the dependence of the warping amplitude on the Fermi energy manifests itself in 
a strong dependence of the diffusion constant on this Fermi energy, leading to several 
important experimental consequences. Moreover, the combination of a strong warping 
with an in plane Zeeman effect leads to an attenuation of conductance fluctuations in 
contrast to the situation of unwarpcd Dirac surface states. 
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1. Introduction 

Topological insulators (TIs) constitute a new state of matter with an insulating bulk 
and an odd number of metallic Dirac cones at their surface [HE]- It was predicted 
[3] and confirmed by angle-resolved photoemission spectroscopy (ARPES) jU El E] that 
the compounds Bi 2 Se3, Bi 2 Te3 and Sb2Te 3 possess such a single surface state (SS). 
Owing to this spin locking property, SSs of TIs are unique metallic systems for the 
study of fundamental magnetotransport properties and for the realization of future 
spintronics based devices. Unfortunately such transport experiments have faced a major 
difficulty so far. Indeed most TIs samples present enough residual bulk conductance 
to overwhelm the actual surface contribution [TJ [8J [9J [ID] . The Shubnikov-de Haas 
oscillations reported in Bi 2 Se 3 crystals originate exclusively from the bulk 3D bands 
indicating a low mobility for the surface Dirac fermions [TJ El E] . Nevertheless recent 
progress has been achieved by using ultrathin films of Bi 2 Se 3 [HI [121 H3] or Bi 2 Te 3 [H] 
exfoliated on high-k dielectric insulator in order to improve the surface/bulk conductance 
ratio and to allow the gating of the surface state. From the gate voltage dependence of 
the total conductance, it was possible to separate bulk and surface contributions using 
a classical two-carrier model [TU [T2] • Most recently the surface conductance of strained 
HgTe samples has been reported in truly 3D slabs with thicknesses exceeding 100 nm 
and in the absence of significant bulk conductance [15]. At low temperatures and for 
small samples, the quantum correction to transport was recently probed experimentally 
: the dependence of conductance on weak magnetic fields displays the expected weak 
antilocalization while the Universal Conductance Fluctuations (UCF) were observed 
in some of those experiments [HJ [TTJ [TS]. Interestingly, the crossover between the 
symplectic and unitary universality classes upon breaking time-reversal symmetry has 
been observed in ultrathin samples of Bi 2 Te 3 [T6j . 

Theoretically, describing the transport properties of topological insulator surface 
states amounts to consider the diffusion properties of two dimensional Dirac fermions. 
Indeed, at the lowest order in k.p theory the SS Fermi surface is circular with a spin 
winding in the plane. At low energy, this two dimensional conductor shares a lot of 
similarities with graphene but with the following important differences: the momentum 
is locked to a real spin as opposed to a A-B sublattice pseudo-spin in graphene, and 
it has a single Dirac cone as opposed to the four-fold degeneracy of the Dirac cone 
in graphene. For this simpler Dirac metal, the conductivity and the induced in-plane 
spin polarization were calculated as functions of the 2D carrier concentration [T7j. Far 
from the Dirac point, surface spin-orbit may generically produce a significant hexagonal 
warping (HW) of the spin texture. As a result the Fermi surface exhibits a snowflake 
or a nearly hexagonal Fermi surface depending on the carrier density, and the spin 
gets tilted out of the plane [18]. Those effects have been confirmed by ARPES and 
scanning tunneling microscopy (STM) experiments performed on Bi 2 Te 3 crystals where 
HW is particularly strong PU H2 HO E21 E3] , and also in Bi 2 Se 3 [21 125]. Characterizing 
quantitatively the HW and therefore the amount of out-of plane spin polarization is a 
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crucial issue in view of potential spintronics applications, and has been the subject of 
ab-initio studies [26] and ARPES experiments [22J . This hexagonal warping is specific to 
the Dirac metal at the surface of topological insulators, and different from the trigonal 
warping encountered at high energies in graphene. Its amplitude is stronger in the 
metallic regime of high Fermi energy which is of interest in the present paper : a full 
description of transport properties in this regime necessitates to include its presence in 
the model. It has been shown fI7\ fIE\ that hexagonal warping enhances perturbatively 
the classical Drude conductivity with respect to the unwarped model [T7j . However, 
as we will show in this article a description of the effects of this hexagonal warping 
requires to go beyond the previous perturbative descriptions. In this paper, we develop 
a description of the diffusion of these surface states which is non-perturbative in the 
warping amplitude. 

We investigate theoretically the classical 2D charge diffusion within an hexagonally 
warped surface states using both a standard diagrammatic formalism and a Boltzmann 
equation approach. In this classical regime, we find that the hexagonal warping on 
one hand strongly reduces the density of states and on the other hand increases the 
diffusion coefficient for the Dirac surface states. The combination of both effects is found 
to correspond to an increase of the classical conductivity as a function of the warping 
amplitude. Interestingly, as the warping amplitude increases with Fermi energy, that 
dependence on the warping amplitude manifests itself in the behavior of the conductivity 
as a function of the Fermi energy. In the second part of this article, we focus on the 
quantum correction to the conductivity, relevant in a phase coherent conductor. These 
quantum corrections are known to depend on the symmetry class of the associated 
Anderson localization problem. Similarly to graphene with only intra-valley disorder, 
the single Dirac cone model is known to correspond to the symplectic / All class. Within 
this symplectic class the Dirac model differs from graphene with long-range disorder 
or electrons with spin-orbit randomness by the presence of a topological term in the 
associated field theory [29, 30J. This topological term explains the absence of Anderson 
localization for these topological insulators surface states [3TJ E2J US] , however it plays 
no role in the diffusive metallic regime. In this regime, the universal properties of both 
weak localization and conductance fluctuations are just the same as for any model on 
the two dimensional symplectic / All class. We will recall these results, and show that 
the dependence of the diffusion coefficient on warping manifests itself in a dependence of 
the weak localization correction and conductance fluctuations away from their universal 
values. In particular, this diffusion constant parameterizes the universal cross-over from 
the symplectic / All class to the unitary / A class when a magnetic field is applied 
perpendicular to the surface states. The shape of the associated cross-over functions 
now depends strongly on the warping amplitude, and thus the Fermi energy. Finally, 
we identify a property unique to the topological insulator surface states : a reduction 
of the conductance fluctuations by the application of a in-plane magnetic field. 

The paper is organized as follows. We present the model in section [2] The classical 
Drude conductivity in presence of HW is derived using a Boltzmann approach (section 
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[3]) and a diagrammatic method (section [4]). In section |5j the quantum corrections to 
the average conductivity and the Universal Conductance Fluctuations are respectively 
derived. The results and conclusions relevants for experiments are summarized in the 
conclusive section [6] 

2. Model 

2.1. Hamiltonian and Warping Potential 

We consider the Hamiltonian describing a single species of surface states of a strong 
topological insulator 

H = hv F a-k + ^a z (k 3 + + k 3 _) + V(r), (1) 

where a = (o~ x ,o~ y ) is the two-dimensional vector of Pauli matrices describing the 
physical in-plane spin of the electron, and k = (k x , k y ) is its in-plane momentum operator 
which satisfy the usual commutation relation with the position operator [r a , kp] = iS a p 
(a, (3 = x,y, z). Note that in Eq.Q, we used for sake of simplicity the dispersion relation 
a-k for the Dirac part instead of the usual axk hamiltonian : switching from one to the 
other amounts to perform an in plane 7r/2 rotation around z in the spin space and has 
no consequence on the following discussion as it does not affect the hexagonal warping 
potential. Concerning the description of the disorder potential V(f), we will follow the 
standard description by considering a Gaussian random potential V(f), characterized 
by a zero average value (V(f)) = and variance (V(r)V(r')) = 7<5(r — r'), where (. . .) 
represents the average over different disorder realizations. 

Dispersion Relation 




Figure 1. Fermi surface in presence of hexagonal warping for various Fermi energies 
with parameters relevant for Bi2Te3 : A = 250 eV.A 3 and vf = 2.55 eV.A and Fermi 
energies from Ep = 0.05 eV to 0.3 eV, corresponding to a warping parameter b ranging 
from to b ~ 0.7. 
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The hexagonal warping term if w = %cr x (k+ + k 3 ^), with k± = k x ±ik y , breaks the 
full U(l) rotation symmetry of the Dirac cone down to a trigonal discrete C3 symmetry 
[T8| 134"] . Time reversal symmetry leads to an additional hexagonal C& symmetry of the 
Fermi surface, hence the name hexagonal warping (Figjl]). The dispersion relation in 
the absence of disorder with k = (k x , k y ) = k(cos9,sm9) : 

e 2 (k) = h 2 v 2 F k 2 + X 2 k 6 cos 2 (3#), (2) 

is obtained by squaring the Hamiltonian Eq.Q with V(f) = 0, and leads to the 
snowflake shape of equal energy surfaces at high energies : see Figjl} Defining 
E F = UvpkF and k = kpk(9), the shape of the Fermi surface at energy E F is defined by 
the equation 

1 = k 2 {9)+4b 2 ~k 6 (9) cos 2 (39). (3) 

Hence the warping of this Fermi surface is naturally characterized by the dimensionless 
parameter 

XE ? (4) 



2fr?v F ' 

which will play a crucial role in the remaining of this paper. This parameter is related to 
geometrical warping deformation factor w introduced in ref. |22j, and defined from the 
maximum and minimum momenta k max , k m i n of a constant energy contour (see Figjl]) 
as 

~ w k max k mm 2 -|- Vo . . 

w=- = - — — with w max = - /= ^ 13.9 (5) 

" J max "'max ~r ^min Z — V <J 



With the notations of Eq. O), we have k min /k max = k(9 = 0), leading to the expression 



Vw(l + w) 2 (w max + w) 2 

" {W) = (!-»)' = ^(«W-«f (6) 
Two values appear remarkable : w = 0, b = correspond to a circular Fermi Surface 
; w — 1,6 = 2/(3\/3) ~ 0.38 to a hexagonal Fermi surface; while w > 1 indicates a 
snowflake shape. 

Using the experimental values for the Bi 2 Se3 compound A = 128 eV./l 3 , vf = 
3.55 eV.A from [21] we obtain relatively small values of warping 0.04 < b < 0.09 
for energies 0.05eV < E < 0.15eV, and similarly 0.0 < b < 0.04 with the values 
A = 95 eV.i 3 , v F = 3.0 eV.i and O.OOeV < E < 0.15eV from [25]. On the other hand, 
the experimental values for the Bi 2 Te3 compound A = 250 eV.A 3 and v F = 2.55 eV.A 
[5j [19] lead to a warping factor ranging from b = 0.13 for E = 0.13 eV to b = 0.66 
for E = 0.295 eV. As we will show below, these large values of warping amplitude are 
indeed beyond the reach of a perturbative approach, and require a non-perturbative 
description of diffusion in the presence of warping. 
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3. Classical Drude conductivity: Boltzmann equation 

The momentum relaxation rate in presence of warping can be obtained from a simple 
application of Fermi's Golden Rule: 

-= 1 ^2n\(k'\k)\ 2 1 8(e(k')-E F ), (7) 

where (k'\k) is the overlap between the eigenstates labelled by k and k' at the Fermi 
surface. The disorder strength is characterized by 7 = UiV 2 , and the dispersion e(k) is 
given by Eq. 

At equilibrium, the occupation numbers are given by the Fermi distribution 
rzi?(e(/c)). In presence of a finite external electric field E, the occupation number function 
f(k) has to be determined by solving the Boltzmann equation: 

- e£.^P = / ^2n ni \(k'\V\k}\ 2 5(e(k') - e(k)) \f(k>) - f(k) 
dk J V"k) 2 L 

This equation is nontrivial to solve due to both the warped shape of the Fermi surface 
and the anisotropy of the scattering. Indeed, even if the scalar disorder is isotropic and 
spin diagonal, the spin-locking property makes the scattering anisotropic (even for a 
circular unwarped Fermi surface). This anisotropy is further increased by the hexagonal 
warping of the Fermi surface. This is apparent when considering the scattering 
amplitude which appears in Eq.|8| : f(9,9' — 9) = \(k\V\k'}\ between eigenstates of 
the disorder- free Hamiltonian H + if w in Eq. ([!]), as a function of the polar angle 9 on 
the incident state \k) and polar angle 9' of the outgoing state \k'). Due to the presence of 
hexagonal warping, this scattering amplitude is highly spin-dependent and anisotropic 
along the Fermi surface, especially at large Fermi energy where the warping factor b is 
large. Moreover this amplitude f(9, 9' — 9) becomes dependent on the incident state 
direction 9, and not only on the relative angle, (9' — 9) as opposed to the case of Dirac 
fermions without warping. These properties are represented in Fig. [2] 
To solve the equation (JsT) , we use the linearized ansatz: 



f(k)=n F {e(k)) + d ^f(e), (9) 
where the angular function f(9) is linear in electric field. We find: 

- eB % = 2i °7 s? m?k " W m - m] • (io) 

where the Fermi wavevector kpgi and the derivative {de/dk)gi are both evaluated at the 
point of the Fermi surface labelled by the angle 9'. The solution of this integral equation 
at fourth order in b is: 

f(6) =ev F T e E x \2cos9 + b 2 (l8cos9 + 5cos(59) - cos (7 9)) (11) 

b 4 1 
- — (522cos# + 20cos(5#) +44cos(7#) +49cos(ll#) - llcos(130)) . 
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Figure 2. Scattering amplitude |(fc|T^|fc')| 2 plotted as a function of the relative angle 
9' — 9 between the incident state k of polar angle 9 and the scattered state k! of polar 
angle 9' . Three different incident directions are plotted to demonstrate the dependence 
on the polar angle of the incident state : 9 — 0, 7r/12, tt/6, while the results are invariant 
by 7r/3 rotation of the incident state. The results are shown for increasing energies 
warping factor b = — 0.7 corresponding to increasing Fermi energies in a given 
material. 



The current is then calculated as: 



jx = I J^^JiW^) - E F ) = aE x , (12) 



leading to the conductivity: 

a = a (Q ) (1 + 8b 2 - 58b 4 + o(6 4 )) . (13) 

The conductivity and density of states (at Fermi level) in the absence of warping (b = 0) 
are given respectively by cr^ = e 2 p^VpTe°^ and = 27T fj v 2 ■ Similarly, one can also 
derive the perturbative (and non-perturbative) expression for the density of state p : this 
is done in Appendix A| and the result is represented on Fig. |3j However, to go beyond 
the above perturbative expansion for the conductivity, a diagrammatic approach turns 
out to be more convenient. Hence we proceed below by developing such an approach, 
comparing its results with the perturbative expansion of Eq. (13). 



4. Classical Drude conductivity : diagrammatic approach 

In this part, we investigate within the standard diagrammatic framework the diffusive 
transport of Dirac surface state with an arbitrary large warping deformation (see 
Eq. ([T])). We start by calculating the single particle Green function averaged over 
disorder whose imaginary part yields the density of states. Then we evaluate the classical 
Drude conductivity from the bubble diagram containing two dressed current operators 
linked by two disorder averaged Green functions. We discuss the dependence of the 
classical conductance as a function of the 2D carrier density and the HW coupling 
strength. 
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Figure 3. Renormalization of the density of states p(Ep) by the warping of the 
surface states, non perturbatively in the warping amplitude b. The resulting correction 
is represented p(Ep)/p^ where p^ = p{\ = 0) is the density for Dirac Fermions 
without warping. 



4-1. Averaged Green's function and elastic scattering time. 

The averaged Green's function is obtained by calculating the self-energy correction and 
averaging over disorder [35J. The corresponding disorder averaged propagator reads: 

(E ± ih/2T e )± + hv F k . a + Xk 3 cos(39)a z 



(G R/A (k)) 



(14) 



(E ± ih/2r e ) 2 - h 2 v 2 F k 2 - X 2 k 6 cos 2 (39) 
where the elastic scattering time for the warped (resp. unwarped) conical Dirac fermion 
r e (resp. r e (0) ) is defined by : 

h 



(0U0) 



pr e = P x 't ( 



7T7 



(15) 



We assume that the phase coherence time ^> r e in order to be in the regime of coherent 
transport. By using the parameterization k(9) of the Fermi surface shape introduced in 
Eq. (|3]), we obtain the non perturbative expression for the the density of states : 



,(o) 



a(b) 



2tt 



dO 



T, 



(0) 



(16) 



P w Jo 2tt 1 + 12 b 2 k\9) cos 2 (39) 

We find an increase of the scattering time r e by the warping. Hence in the presence 
of hexagonal warping, the surface states of topological insulators are relaxing carrier 
momentum less effectively than without warping. This is a manifestation of the 
increased anisotropy of the scattering amplitude by the hexagonal warping correction, 
illustrated in Fig|2| Indeed in the presence of hexagonal warping, the neighboring states 
of an incident Dirac fermions have smaller overlaps than states that are ±27r/3 apart, 
which effectively induces a very anisotropic scattering. Nevertheless from this result, one 
can not yet conclude that the Drude conductance is enhanced as the density of states at 
the Fermi energy is also renormalized by the warping. Indeed the Drude conductance is 
determined by the diffusion constant D and the transport relaxation time rather than 
the elastic scattering time r e governing momentum relaxation. To access this classical 
conductivity, we follow the standard diagrammatic procedure in the following section. 
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4-2. Diagrams for the classical conductivity 

For a given impurity configuration, the conductivity is given by the Kubo formula : 

h 



2ttL ; 



. )MV [j x G R j x G A ] 



(17) 



where we have neglected subdominant contributions involving products G R G R and 
G A G A %5jg\ • In this expression and below, Tr denotes a trace over the electron's Hilbert 
space (momentum and spin quantum numbers). The current operator j a is obtained 
from the Hamiltonian Eq. by inserting the vector potential via the minimal coupling 
substitution k — >■ k — {e/K)A and using the definition j a = j^-: 



3A 



;is) 



Note that this current contains the usual — ev F a x term for linearly dispersing Dirac 
system, and an additional term quadratic in momentum originating from the HW. 

After averaging over all impurity configurations, the classical mean conductivity 
is obtained from the bubble diagram of Fig. [4] where the propagating lines represent 



the retarded and advanced disorder averaged Green functions Eq. (14). This classical 



conductivity is the sum of two terms : a a and o~b- The contribution of diagram A in 



k + f , E F + hcu 




+ 



v - ' 

(A) 




Figure 4. Diagrammatic representation of the two contributions a a and erg to the 
classical conductance 



Fig. |4] is defined at q = and u> = by 



h 



0~A 



dk 



tr 



Jx (k)(G R (k,E F ))j x (k)(G A (k,E F )) 



(19) 



2tt J (2tt) 2 

where tr denotes a trace only over the spin quantum numbers. Performing explicitly 
the trace and the integration over the momentum k = \k\, we obtain the following 
expression: 

he 2 v 2 F a{b) + 5/3(6) + 5(b) 



OA 



27T7 



a(b) 



(20) 



which is non perturbative in b. The function a(b) has been defined in Eq.(16) and we 
have introduced : 

" 27T d6 1 



a(b) 



2tt1 + 12 b 2 k\6) cos 2 (3#) 



(21) 
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Figure 5. Schematic representation of the Diffuson and Cooperon structure factors. 
The dotted lines correspond to disorder correlators, while the plain and dashed lines 
represent retarded and advanced Green's functions. 

f**dB 4fe 2 cos 2 (30)fc 6 (0) 
Jo 27U + 126 2 cos 2 (30)^(0)' 

pdQ 36b 2 (fe 4 (0) - fe 6 (0)) 
Jo 27rl + 126 2 cos 2 (30)^ 4 (0)' 

where /c = fcpfc(0) with fc(0) introduced in Eq. (J3J) - 

The contribution of the diagram B in Fig. [4] accounts for the contribution of the 
so-called Diffuson [36J. Retaining only the dominant contribution of the integral, we 
can write it explicitly as 

a B {q,u) = A tr [j x T^(q,u)J x ] , (24) 
where the vertex operator J is defined by 

/dk - 
— (G A (k,E)) Jx (k) (G R (k, E)) = J o x . (25) 

A contraction over the spin indices is assumed, resulting in the proportionality to o~ x . 



The Diffuson structure factor in Eq.(24) is defined diagrammatically in Fig. [5J 
and satisfies the usual recursive Dyson equation, solved by the expression F^(q, u) = 
7 [l ® 1 — 7P^(g, a;)] where is the polarizability 

/dk 
j^y 2 (G R (k, E)) g) (G A (k -q,E- co)}, (26) 

where £g> denotes a tensor product of the spin Hilbert spaces. The resulting structure 
factor is naturally decomposed into 4 spin modes. However, as opposed to the 
case of non-relativistic electrons [US], for finite but small q, the non-diagonal character 
of the Dirac Green's functions leads to unusual terms in this structure factor such 
as (q.cr) £g> (q.a), (<f.<?) <8> 1. We recover the standard singlet and triplets states only 
in the diffusive q — y limit. In this limit, the only gapless mode is the singlet 
state, characteristic of the symplectic / All class. While this mode and the associated 
Cooperon singlet determines the quantum corrections to diffusion, it does not contribute 



to the classical conductivity of Eq. (24). The only classical contribution comes from one 
of the massive triplet states. In full generality, the q — > limit of the structure factor 
can be parameterized according to 

= ail ® 1 + a 2 (a x <g> a x + a y g> a v ) + a 3 a z <g> a z . (27) 
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Using the parameterization of Eq. (25) and performing the resulting trace we obtain 

He 2 v 2 F (a(b) + 2/3(b)) 2 



°B = J (Ol - 0,3) 



2tt7 a(b)(a(b) + (3(b)) 



(28) 



Adding the two contributions of Eq. (20) and (28) leads to the full expression for the 
classical conductivity : 

o~a + o~ B = e 2 pD. 



(29) 



The evolution of this classical conductivity as a function of the warping amplitude is 
represented in Fig. |6l Moreover, this expression of the conductivity allows us to define 



cr 



cr(b = 0) 



./'order b 2 




0.1 0.2 0.3 0.4 0.5 0.6 



Figure 6. Longitudinal conductivity as a function of the parameter b. Solid line is 
the exact calculation, dotted line is the 4th order expansion in b calculated through 
the Boltzmann equation, dashed line is the 2nd order. 



the diffusion constant D from the Einstein relation, whose final expression reads 



D 



a(b) + 5/3(6) + 5(b) (a(b) + 2/3(b)) 2 \ 
a(b) + a(b)(a(b) + (3(b)) ) 



(30) 



which is non perturbative in the warping parameter b. In the limit b — > of absence 
of warping, we recover the known result for Dirac fermions in the presence of a scalar 
disorder with a diffusion constant D = v 2 F r e = VpTtr/2 and a transport time r tr = 2r e 
accounting for the inherent anisotropic scattering of Dirac fermions on scalar disorder. 
The renormalized diffusion constant is found to increase as a function of b (Fig. [7]). 
This dramatic increase of D within the range of experimentally relevant warping 
parameter b signals an effective strong increase of the anisotropy of scattering when 
hexagonal warping is present. Note that in the present case, this diffusion coefficient 
accounts for both the renormalization of the averaged Fermi velocity and the scattering 
amplitude corresponding to a renormalized transport time. Moreover, the comparison 



of the expression of Eq. ( 30 ) with its perturbative expansion to order b A demonstrates 



that the experimental values of warping b are beyond the reach of this perturbative 




Figure 7. Evolution of the diffusion constant as a function of warping intensity b. The 
solid line corresponds to the expression non perturbative in the warping amplitude b, 
while the dotted and dashed curves correspond to perturbative results to order b 2 and 
b 4 . 



expansion (and the similar perturbative studies of [27J[28]) : to describe accurately the 



experimental situation at high doping, the full exact expression Eq. (30) is required. 



5. Quantum corrections to the conductivity 

In this section, we focus on the regime of quantum transport corresponding to a phase 
coherent diffusion of the surface states which can be reached for small sample size and at 
low temperatures. This regime corresponds to the situation where the Fermi momentum 
kp satisfies the condition kpl e S> 1, which naturally corresponds to the situation where 
the warping of the Fermi surface is strong. We focus on the first two cumulants of the 
distribution of conductivity and compute diagrammatically the corresponding quantum 
corrections non perturbatively in the warping amplitude b, extending the results of Refs. 
[271 E8]. 

5.1. Universality classes 

The study of non-interacting metals perturbed by disorder corresponds to the well- 
studied problem of Anderson localization of electronic waves. In this framework, the 
transport properties depend on both dimension and universality classes determined by 
the symmetries preserved by disorder. These universality classes encodes both the 
universal properties of the Anderson transition, but also the universal properties of 
the weak disorder metallic regime of interest experimentally in the present case. The 
model considered in this paper H = Uvfo ■ k + V (r) is described at large distances by 
the standard All / symplectic class : time-reversal symmetry T is preserved, but due to 
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the momentum-spin locking it squares to T 2 = —Id [371 [38] . This is also the Anderson 
universality class describing non-relativistic electrons in the presence of random spin 
orbit disorder. However while both models possess the same universal weak localization 
properties, they differ in their strong disorder behavior : the topological nature of the 
single Dirac species in the present model manifest itself in the appearance of a topological 
term in the field-theoretic description of the d = 2 All class [291 130] . The presence of 
this topological Z 2 term modifies the strong disorder / low conductance behavior of 
the All class [311 E2J [33] , but plays not role in the weak disorder regime. Hence the 
behaviors of the first two cumulants in the weak localization regime of the single Dirac 
model with scalar disorder are exactly those of the standard All / symplectic class. 
Moreover adding the HW term H w does not break Time Reversal Symmetry, and thus 
does not change the universal symmetry class. In contrast, either a magnetic disorder 
or a magnetic field introduced through a Zeeman or orbital term in the Hamiltonian 
breaks time-reversal symmetry and induces a crossover from the All/symplectic to the 
A/unitary class characterized by the absence of time reversal symmetry. Such a crossover 
between the symplectic and unitary ensembles has already been observed by depositing 
Fe impurities at the surface of ultrathin samples of Bi 2 Te3 [32]. In this experiment, 
the exchange field of the Fe atoms induces the Zeeman coupling for the SS carriers and 
destroys the weak antilocalization signature. Nevertheless the extreme thickness of the 
sample is likely to yield an important hybridization between the bulk, and top/bottom 
SSs while here we consider the response of a single isolated SS. 

One of the characterizations of Anderson universality symmetry classes is the 
number of independent low energy modes in the diffusive regime. Those modes 
correspond to the various Cooperon/Diffuson modes which indeed parametrize the target 
space of the underlying non-linear sigma field theory. The symplectic class corresponds 
to two modes, one singlet Diffuson mode and one singlet Cooperon mode [SB]- The 
breaking of time reversal symmetry by e.g. an orbital magnetic field will suppress the 
Cooperons whose propagator acquires a mass and is no longer diffusive : this induces 
a cross-over from the symplectic class to the unitary class, characterized by a single 
diffusive mode. This cross-over occurs on a length-scale determined by the mass acquired 
by the Cooperon singlet. The suppression of the Cooperon is attributed to destructive 
interference between time-reversed paths, accounted for by the retarded and advanced 
Green's functions in the Kubo formula. Such a cross-over will be discussed in section 

El 

5.2. Weak antilocalization 

When winding around the Fermi surface, the spinorial electron wave function acquires 
a 7r Berry phase which is responsible for weak antilocalization (WAL) phenomenon. 
In graphene, the presence of intervalley scattering leads to a crossover from weak 
antilocalization to weak localization when the ratio of intervalley over intravalley 
scattering rates is increased [35] 139] . In contrast, TIs with a single Dirac cone provide 
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ideal systems to measure WAL once the issue of the spurious bulk conductance is solved. 
Actually some experiments have already reported a strong WAL signal in thin films and 
in 3D HgTe slabs [15] . 

The quantum correction to the conductivity in the class All corresponds to a weak 
anti localization correction : its determination within diagrammatic theory for Dirac 
fermions is recalled in Appendix B In the limit of a phase coherent sample of size 



Ls C I where is the phase coherent length, it reads: 



fe 2 \ f 1 



2 



<fa>= UM^ = ^ ln(V4) (31) 

where i e is the elastic mean free path. This result is independent of the particular 
model within the class AIL However, it depends on the the amplitude b of warping at 
high Fermi energy through the diffusion coefficient entering the phase coherence length 
L<f, = a/-D(&)t0, where is the phase coherence time. 

When a transverse magnetic field is applied, the standard derivation to account 
for the magnetic orbital effect of this result still holds |36J: it amounts to derive the 
Cooperon contribution from the probability of return to the origin of a diffusive path in 
the presence of the magnetic field. The corresponding correction to the conductivity is 
thus given by [36] : 



(5a(B)) 



e 2 



2 B 12 B 



(32) 



where the characteristic fields B e = h~/ (4eD(b)r e ) and B^ = h~/ (AeDib)^) have been 
introduced, and where \1/ is the Digamma function. The diffusion constant D(b), non 



perturbative in b, is given by Eq.(30). This corresponds to the result obtained for 
graphene when inter- valley scattering can be neglected (35] . In our case, this expression 
implies that in a given sample, the shape of the weak anti-localization correction as a 
function of B will evolve as the Fermi energy is varied within the Topological Insulator 
gap and the warping amplitude b is varied. This effect will be discussed further in 
section |6l 

5.3. Universal Conductance Fluctuations 

As explained above, in the absence of magnetic field the Universal Conductance 
Fluctuations (UCF) result for the symplectic class results still holds and one finds, 
similarly to graphene [10] : 



< faJ > = 12 (9 vf t 7- (33) 

The derivation of this result for Dirac fermions is recalled in |Appendix C| Defining the 
phase coherent = 'D(6)r^(T) and thermal lengthscales L T = ^h~D(b)/T, we can 
focus on different regimes : for L -C L<a, Lt, a proper regularization of the integral in 



(33) leads to the universal value 



M -?(i)£S5R? aUM (i)' (34) 
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This result is independent of the diffusion coefficient D(b) and thus independent of the 
warping amplitude. On the other hand, in the other limits [36] : 



(&t 2 )~-(-1 U^l for^«L,L T , (35) 
MVe 2 \ 2 , / L 



, j y h j In ^— ) for L T < L < L^, (36) 

rX 2/2\2 / r 



? J h J In IjZJ for L T < L < L. (37) 

Hence a strong dependence of these UCF on the Fermi energy through the warping 
amplitude b is found in all these cases. The introduction of a transverse magnetic field 
induces a crossover from the symplectic class to the unitary class, where the amplitude 
of the fluctuations is reduced by a factor of two. This crossover is described as 



(6a(Bf) = l -(5a 2 ) 



l + ^'(- + — 
B v 2 B 



(3* 



where \I/ is the Digamma function and B^(b) = h/AeD(b)r^. The dependence in b of the 
diffusion constant will affect the value of : this characteristic field for the suppression 
of the WAL correction and this reduction by a factor 2 of the conductance fluctuations 
will decrease when the Fermi level is raised away from the Dirac point. 

5.4- In-plane magnetic field: interplay between Zeeman and warping effects 

We now consider the effect of an in plane Zeeman field on the transport properties of 
these Dirac fermions. This can be accounted for by adding to the Hamiltonian ffl a 
term 

H z = 9^b {o- x B x + o-yBy) . (39) 

Without warping, such a Zeeman field acts exactly like a constant vector potential and 
can be gauged away. This amounts to shift uniformly the momenta by —gjj, B B /(hvp). 
Hence without warping, an in-plane Zeeman field does not modify the scattering 
properties of an (infinite) disordered sample. However the presence of the hexagonal 
warping term in (fl]) breaks this invariance : the shape of the Fermi surface is now 
modified by the Zeeman field. As a consequence, the scattering amplitudes f(9, 6' — 6) 
acquire also a B dependence. We can naturally expect that this Zeeman field which 
modifies the scattering amplitudes redistributes the scattering matrix of the samples, 
and leads to conductance fluctuations induced by an in-plane magnetic field. 

To describe quantitatively the effect of this Zeeman field, we extend the above 
diagrammatic analysis perturbatively in B = g^sB/Ep . This Zeeman field, which 
breaks time reversal symmetry, induces a crossover from the symplectic to the unitary 
class. This crossover can be accounted for by the evolution of the Cooperon structure 
factor : its singular part does not correspond anymore to a diffusive singlet component 
but to the diffusion of a massive singlet: 

T c (Q) = -—^— 1 -\S}(S\. (40) 
Dq z r e + m(o, B) 
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The mass m(b, B) which encodes the effects of the Zeeman field, is calculated 
perturbatively to the second order in B : m(b, B) = m(b)B 2 . The non-perturbative 
result is represented in Fig. ^ We recover the expected result m(0, B) = m(b, 0) = 
corresponding to the single singlet diffusive Cooperon mode of the symplectic class. This 




Figure 8. Evolution of the massive term m(b) in the diffusion equation of the 
Cooperon, in presence of the hexagonal warping and in-plane Zeeman magnetic field, 
as a function of the warping parameter b. 



mass is associated with a length parameterizing the cross-over from the symplectic to 
the unitary class, defined by Lb = \J DT e /m(b, B). Beyond this length scale, we recover 
a standard unitary weak localization, and conductance fluctuations are reduced by a 
factor 2 : 

1 



(5a(B)) 



(Sa 2 (B)) = 6 



1 



e 
h 



1 

V 



f(L/£ e ,L/L) 



(41) 



(6a 2 (B 



(<f + ^ 2 ) 2 

6)) + f 2 {L/L), 



+ 



(q 2 + L B 2 + L 



-2\o 



(42) 



where L is the longitudinal size of the topological insulator surface, is the phase 
coherence length and L is defined as L~ 2 = L^ 2 + L^ 2 . The function / and f 2 depends 
on the geometry of the sample. In the case of an infinite 2D sample, one recovers only 
a logarithmic dependence on Lji e for /. For a finite size sample, one has to to replace 
the integral by a sum over the compatible wavevectors [36j . 

As expected, we find that the UCF extrapolates between the value for the 
symplectic and the unitary class, over the length scale LB(b,B). An estimation of 
the evaluation of this length as a function of the magnetic field or the warping term is 
given by : 

c(b) 



L 



B 



'D{b)/D(b = 0)m(6) 
B 



B 



(43) 
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where the evolution c(b) with the warping amplitude is shown in Fig. [9] For experimental 
realistic values, b ~ 0.7, c(0.7) ~ 1, Ep ~ 0.3 eV and p// B ~ 5.10" 4 eV.T" 1 . This 
gives for a magnetic field around 1 T a characteristic length to observe the crossover 
Lb — 1000 £ e which is is rather large but within experimental reach. 

c(b) 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 9. Dependence in b of the characteristic length to observe the crossover from 
symplectic to unitary class, for an in-plane Zeeman magnetic field. 



6. Summary of Results and Conclusions 

In this section, we summarize the main results of our approach focusing on the 
experimentally relevant aspects. In particular, we emphasize the effect of strong 
variation of the warping amplitude b as the Fermi energy is varied in a given material. 
Universal properties of transport in the quantum regime are not affected by the presence 
of this warping. However, both the departure from these universal properties as well 
as the classical regime are entirely characterized by the diffusion constant D(b) which 
itself depends in the warping amplitude. 

The transport properties in the incoherent classical regime are characterized by both 
the density of states p{Ep) and the diffusion constant D. The dependence on the Fermi 
energy of both quantities is strongly affected by the presence of the warping of the Fermi 
surface, whose amplitude b itself depends on E F . This is shown on Fig. [TUJTTJfor two sets 



of parameters corresponding to the Bi 2 Se 3 and Bi 2 Te 3 compounds. As a consequence of 
these results, the conductivity a acquires a strong Fermi energy dependence shown on 
Fig. [TTJ which can be directly probed experimentally. 

Similarly, the quantum corrections to transport depends strongly on the Fermi 
energy through the dependence of the Diffusion coefficient D on the warping amplitude. 
Indeed, the shape of the typical measurement of the weak (anti-) localization correction 
through the dependence of the conductivity on a magnetic field perpendicular to the 
surface depends solely on the diffusion constant. This diffusion constant depending on 
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D (A = 0) 




Figure 10. Dependence of the diffusion constant on the Fermi energy for the values 
of A, vf corresponding to the Bi 2 Se3 and Bi 2 Te3 compounds. The results depend on 
the amplitude of disorder parameterized by the mean free path l e . Here, they are 
normalized with respect to the diffusion constant in absence of warping D{\ = 0) 
which is independent of the Fermi energy and incorporates the dependence on l e . 




£ F (eV) 



Figure 11. Dependence of the density of states (Left) and classical conductivity 
(Right) on the Fermi energy for values of A = 128 eV.A 3 , vp — 3.55 eV.A for Bi 2 Se3 
and A = 250 eV.A 3 and vp = 2.55 eV.A corresponding to the Bi 2 Te3 compounds. 
On the left figure, the dotted lines corresponds to the standard linear density of 
states for Dirac fermions without warping A = 0. On the right figure, the results 
are represented as a ratio with the conductivity in the absence of warping <r(A = 0) 
which is independent on the energy and incorporates the dependence on the disorder 
strength. 



the Fermi energy, the associated anti-localization curve depends itself on this Fermi 



energy as shown on Fig. 12 While the amplitude of conductance fluctuations are 
universal in the limit of an entirely coherent conductor, their amplitude in a realistic 
situation where ~ L will be parameterized by a universal function of the diffusion 
coefficient. We have shown moreover that these fluctuations depend in a remarkable 
way on an in-plane Zeeman magnetic field. The amplitude of this effect depends on the 
ratio between the associated magnetic dephasing length Lb and the mean free path £ e . 
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(&r(B)) (&r(B)) 
<<5cr(B = 0)> <<Scr(B = 0)> 

1.0 1.0 




B 



-1.0 -0.5 0.0 0.5 1.0 B„ -1.0 -0.5 0.0 0.5 ~ 1.0 Bo 



Figure 12. Dependence of the weak localization correction (5a(B)) / (5a(B = 0)) 
on the Fermi energy for the values of X,vf corresponding to the Bi 2 Te3 (Left) and 
Bi 2 Se 3 (Right) compounds. We have chosen to scale the magnetic field as B/B where 
Boll = 4>q = h/e to avoid any energy (or warping) dependence of this rescaling field. 
The results show a clear dependence on energy of the magnetic field characteristic of 
weak localization decay. 



The dependence of this ratio on the Fermi energy is shown on Fij 
sets of values used above. 



13 for the two same 




500- 



0.10 



0.15 



0.20 



0.25 



1130^ (eV) 



Figure 13. Dependence of the magnetic dephasing length Lb on the Fermi energy 
for parameters corresponding to Bi2Te3 (Red) and Bi2Se3 compounds. 



In conclusion, we have shown that it is essential to take into account the hexagonal 
deformation of the Dirac cone occurring at the surface of topological insulators such as 
Bi 2 Se3 and Bi 2 Te3 to accurately describe both their classical and quantum transport 
properties. In particular, we provide a formula describing the evolution of diffusion 
constant D for arbitrary strength of the warping. Since this warping amplitude increases 
with the Fermi energy of the surface states, we predict a dependance on doping of 
transport properties different from those predicted for a perfect Dirac cone. 

This work was supported by the ANR under the project 2010-BLANC-041902 
(ISOTOP). JC also acknowledges support from EU/FP7 under contract TEMSSOC. 
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Appendix A. Non Perturbative Density of states 



In the presence of warping, the density of states at the Fermi level p{ep) is modified. 
One has: 

" +o ° kdk r 27r de 



2tt 



-5 



2tt 



a/ (hvk) 2 + \ 2 k & cos 2 (3fl) 



(A.I) 



Performing the angular integration first, the density of states is obtained as the integral: 



dy 



2{7ihvf J y _ v /(i- y )(46V + y-i)' 



where y_ is the unique real solution of the equation 4b 2 y 3 + y 
formula for solving equations of the third degree, we have: 

1 

V- 



(A.2) 

0. By the usual 



3V3b + Vl + 27b 2 - V Vl + 27b 2 - 3V3b 



(A.3) 



2^3b _ 

It is convenient to introduce the parameterization sinh<p = 3\/3b, giving y_ = 
3sinh(<p/3)/ sinhyj. With the change of variables y_ = 1 — 1/t, the integral in (A.2) can 
be rewritten: 



2V2b(irhv) 2 



dt 



y/(t - l) 3 - t 2 /{2b 2 ) 



(A.4) 



Using Eqs. (17.4.70)— (17.4.72) from [H], the integral in Eq. (A.4) can be expressed in 
terms of a complete elliptic integral of the first kind, giving the density of states: 



2K 



l+f sinh 2 (y/3) 
(l+lsinh 2 ^)) 2 ^ V 1 + 4sinh2 ( < ^/3) 



2(nhv) 2 (! + 4 sinh 2(^/ 3 )) 3 / 4 (! + 4sinh 2 (^/3)) 



1/4 



(A.5) 



In the limit of A — >■ 0, cp — > 0, and (A.5) reduces to the density of states in the absence 
of warping. For large A, we have b ~ 2e </p 3v / 3, and the density of states behaves as: 
p(e) ~ 3- 1/i n~ 2 K(l/2 - v / 3/4)(eA 2 )^ 1/3 . Since the density of states goes to zero for 
large and small energy, it has a maximum at a finite value of e. By simple scaling, the 
maximum is obtained for an energy e* = C iy /(hv F ) 3 /X and p(e*) = C 2 (hv F \y 1/2 . 



Appendix B. Weak anti-localization correction for Dirac fermions 

In this appendix, we derive the quantum correction to conductivity for two dimensional 
Dirac fermions in the absence of warping (b = 0). In this case of linearly dispersing 
Dirac fermions, the current operator j x = —evpO x is no longer a function of k. The 
renormalization of this current operator by vertex corrections can be written as 

s> '-' J.,- ■ J,P"Y P 2/,. (B.l) 

where S x stands for the renormalized operator. The quantum corrections to conductivity 
are associated with contributions between interferences of loops of diffusive paths. 
This interference between a loop and its time-reversed contribution is described as the 
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propagation of the so-called Cooperon. Its propagator is defined through the Dyson 
equation T C (Q, u) = 7(1 <g) 1 - 7-P c (<3, w)] where P c is : 



dk 



(G R (k,E)) (G A (Q-k,E-u)). 



(B.2) 



By a time- reversal operation on the advanced component, we can relate the Cooperon to 
the previously identified Diffuson propagator. Only the diffusive modes of the Cooperon 
contribute to the dominant quantum correction, leading to a structure factor for the 
Cooperon : 



T~e DQ 2 

7 1 1 



[1 <g> 1 - a x ® a x - a y ® o y - a z ® o z \ , 



r e DQ 2 4 

where D is the diffusion constant, D = v 2 F r e in the absence of warping. 



(B.3) 
(B.4) 




r c (Q) 



Figure Bl. Diagrammatic representation of the quantum correction to conductivity 
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Figure B2. Diagrammatic representation of the dressing of the Hikami box 



The weak anti-localization correction is obtained by the contraction of a Cooperon 



propagator and a Hikami box, as represented diagrammatically on Fig. Bl This Hikami 



box is the sum of three different contributions represented in Fig. |B2[ We express the 
first of these contributions as 

h 



(5o- ) 



2tt 



-Tr 



G (k)T, x G (k)T (Q)G (Q - k)Z x G A (Q - k) 



(B.5) 



where we used the notations introduced in the article Tr for the trace over all the 
quantum numbers (spin and momenta), tr for the trace over the spin indices, and 
for the trace over the momentum k. Special care has to be devoted to the order 
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of the spin indices in these expressions. In the expression Eq.(B.5), the Q integral 



is dominated by the small Q contribution originating from the diffusive modes of the 



Cooperon. This justifies a posteriori the projection on the single diffusive mode in (B.3 ). 



Focusing on the most dominant part of this expression, we can set Q — > except in the 
Cooperon propagator, the Green's functions being regular in k. Similarly, we can pull 
the renormalized vertices out of the integral and focus on the Hikami box : 



G A {k)a x G R {k) <g> G R {~k)a x G A {-k) 



<k 

p(E F 



h 



16 



[3a x <g> a x + a y <g> a y - 4 1 <g> 1] . 



(B.6) 



Similarly, the two remaining diagrams contributing to the Hikami box in Fig. B2 

expressed as 



arc 



H? = 4 7 




k </tfi 



G A {k)a x G R {k)G H {-q l ) ® G H {-k)G R {q 1 )a x G A {q l ) 



iR/ 



iR( 



71 

16 



p{Ei 



2Te 

h 



[t®t-a x ® a x ] 



(B.7) 



and : 



H: 



c 



47 




k Jqi 



G A {-q l )G A {k)a x G R {k) ® G H {q l )o- x G A {q l )G A {-k) 



1R1 



7T 



p{Ei 



2r P 



[1<S)1- a x ® a x ] 



16'" v * ' v h 

Summing these three contributions we obtain : 



(B.8) 



H 



c 



71 



cr x + a v <g> a y - 2 1 ® 1] 



(B.9) 



2Te 

h J 16 

The resulting weak anti localization correction is obtained via the final contraction with 



a Cooperon propagator as shown on Fig.Bl, and leads to the expression: 



(5a) 



hp(E F )(-ev F ) 2 f2r, 



2tt 



h 



71 

16 



tr (a x g) a x + a v <g> - 2 1 ® 1) / T(Q) 



(B.10) 



Appendix C. Universal Conductance Fluctuations for Dirac fermions 

To derive the conductance fluctuations, we need to take into account two kinds of 
diagrams containing either Cooperons or Diffusons. The Hikami box for Cooperons has 
been calculated previously. Proceeding similarly with the Diffuson instead of Cooperon 
structure factor we obtain the Hikami box for Diffusons : 



H 



D 



p(E t 



2r e \ 71 



[2 181 + 0- 



a' 



a" 



(C.l) 



h J 16 

We have already performed an integration over the momentum (arising from the Kubo 
formula) in these expressions for the Hikami boxes. Hence we only need to plug a 
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Figure CI. Diagram for the conductance fluctuations with Cooperons 




Figure C2. Diagram for the conductance fluctuations with Diffusons. 



Diffuson (resp. Cooperon) structure factor between two H D (resp. H c ). Summing 

(C.2) 



these two diagrams (Fig. C2 and Fig. CI ) we obtain : 
(5a{) = 



2\ 2 



1 f 1 



h) W/g 4 ' 

The second part of the conductance fluctuations come from the diagrams represented 
in Fig. C3 that we have not yet considered. They require the determination of two 
additional Hikami boxes (one for Diffusons and one for Cooperons) : 

p{E F ) (^j ^ [1 <g> 1 + a x <g> o*\ 



D 



H 

H c 



a 



The final results after contraction in spin space of these diagrams is 

-2\ 



(C.3) 
(C.4) 

(C.5) 



Summing the two contributions (]C.2[) and ( |C.5[ ), we finally get the result of Eq 33 

1 



(6a 2 ) = (Sal) + = 12 



2\ 2 



h 



(C.6) 
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Figure C3. Diagrams for the second contribution to conductance fluctuations 



Appendix D. Quantum correction for a warped Fermi surface 



Appendix B| and |Appendix C| show explicitly the derivation of the weak antilocalization 
correction and the conductance fluctuations for Dirac fermions. As expected, the 
corresponding results display no dependence in the only relevant parameter to 
characterize diffusion : the diffusion constant. These results are naturally expected 
to hold when taking into account the hexagonal warping term. To explicitly show this 
independence, we determine the value of the quantum correction to conductivity in the 
general case where the Fermi surface possesses the hexagonal deformation. 

The first step is to obtain the new Hikami box H c , the difficulty arising from the 
dependence of the current operator on the momentum j x = e (—vpa x + ^c z {k x — ky))- 
Recalling that E x = j x + j x P D T D , we express the "naked" Hikami box as: 

H£ = f \G A (k)Y l x G R (k)G R (-k)Y l x G A (-k) \ . (D.l) 

Jk L 

We perform this integral using polar coordinates, and integrate the radial part : 
H c_ e 2 vlrl T e fd6 A + BA + (g) A_BA_ 

° " 4 7 n Te (°)y 2tt1 + I2b 2 k(6y cos 2 (3fl) 
A±=l± k{9) (cos 6a x + sin 6a y ) + 2b ~k\6) cos(3#)a z (D.3) 

B = (2+ J a x + Qb k 2 (6)cos(26)o- z (D.4) 



a + fi s 

The two "dressed" Hikami boxes, where an impurity line links two Green's 
functions can be calculated by the method used above, and we can write the result 

2 2 2 

as H c = 4 ^ e H(b). We need to close this Hikami box with a Cooperon structure 

factor r c = 7" D(b)Q 2 1 (•"■ ® — ° X ' ® ° X ~ a%l ® ° V ~ ° Z ® P er f° rmm g the sum on 
the spin space h(b) = tr[iif(6)| (1 <8> 1 — o x ® a x — a y <8> a y — a z <8> a 2 )] we obtain the 
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quantum correction to the conductivity : 

We found numerically the correction -f^P to be constant equal to 1 and independent 
of b. 
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